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Convcntional Bell and Stirling numbers arise naturally in thc normal ordering of 
simple monomials in boson operators. By extending this process we obtain gen- 
eralizations of these combinatorial numbers, defined as coherent state matrix ele- 
ments of arbitrary monomials, as well as the associated Dobinski relations. These 
Bell-type numbers may be considered as power moments and give rise to positive 
measures which allow thc explícit construction of new classes of coherent statcs. 

1. Introduction 

A defining characteristic of coherent states is the resolution of unity prop- 
erty, which is another expression of the classical Stieltjes moment problem. 
For the conventional coherent states of quantum òptics the moments are 
simply n\, in some sense the simplest combinatorial numbers, leading to a 
simple exponential weight function. Not many other solutions to the mo- 
ment problem are known, and it is in a sense not surprising that one should 
seek for solutions in the realm of other combinatorial integers. One such 
set is the Bell and Stirling numbers, which we define explicitly later, and 
which arise naturally in the normal ordering properties of boson operators. 

The Standard boson commutation relation [a, a^] = 1 can be formally 
realised by identifying a = 4- and a) = x, since = 1. In the present 

note we shall use both of the above forms. Integer sequences arise naturally 
when considering the action of (xd/dx) n on f{x) as the following low order 
examples reveal: 




xf(x) 



xf(x) + 3x 2 f"(x) + x 3 f"(x) etc, 



xf(x)+x 2 f"(x) 



(1) 



(3) 



(2) 
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which in general can be written as 

d_ 

dx 



(x ^) n f(x) = S ( n > k)x k (d/dx) k f(x) (4) 
k=i 



or, alternatively 1 

n 

(a^a) n = ^S(n,fc)(a t ) fe a fe . (5) 
k=l 

The Stirling numbers of the second kind S(n,k) appearing in Eqs.(3) and 
(5) have been known for over 250 years 2 . Eq.(5) exemplifies the normal 
ordering problem, that is, finding the form of (ata) n with the powers of 
a on the right. Although explicit expressions for S(n, k) are known 3 , of 
particular interest here are the Bell numbers B{n) given by the sums 

n 

B(n) = Y,S(n,k), n=l,2,... (6) 
fc=i 

with B(0) = 1 by convention (S(n,0) = S n fl)- A closed-form expression 
for B(n) can be found by considcring the action of (xd/dx) n on a function 
f(x) having a Taylor expansion around x—0, i.e. f(x) — Y^k^o c k xk ■ 
Applying Eq.(4) to f(x) gives: 

, oo 

(x-) n f(x) = J2ckk n x k . (7) 

Now specify f(x) = e x in Eq.(7) so that = 1/fc! and deduce from Eq.(4) 
that 

(voÈ^ = E%^ (8) 

which for x = 1 rcduces to 

oo , n n 

(l/e)Efc!=E^ fc ) = B W· ( 9 ) 

fc=0 ' fe=l 

Equations (8) and (9) are the celebrated Dobinski formulas 2,3 ' 4,5 which 
have been the subject of much combinatorial interest. For complcteness, 
we recali the combinatorial defmitions oíB(n) and S(n, k): B(n) counts the 
number of partitions of a set of n distinguishable elements; S(n, k) counts 
the number of partitions of a set of n distinguishable elements into k non- 
empty sets. Eq.(9) represents the integer B(n) as an infinite series, which is 
however not a power series in n. An immediate consequence of Eqs.(8) and 
(9) is that B(n) is the n-th moment of a (singular) probability distribution, 
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consisting of weighted Dirac delta functions located at the positive integers 
(the so-called Dirac comb) : 



poc 

B(n)= x n W(x)dx, n = 0,l,... (10) 
Jo 



where 



W(x) = (l/e)f: 6 \^. (11) 



k=i 

The discrete measure W(x) serves as a weight function for a family of 
orthogonal polynomials Cn\x), the Charlier polynomials 6 . The exponen- 
tial generating function (EGF) of the sequence B{n) can be obtained from 
Eq.(9) as 

e^" 1 >=f;B(„)^. (12) 

k=0 

This equation is related via Eq.(5) to a formula giving the normal ordered 
form 7 ' 8 of e Aata 

e Aat a = ^(gAa+a) = . ^ a{e* -l) . ^3) 

The symbol Ai denotes normal ordering and : A(<r,a) : means expand A 
in a Taylor series and normally order without taking account of the com- 
mutation relation [a, at] — 1. We stress that in the derivation 7 ' 8 of Eq.(13) 
no use has been made of the Stirling and Bell numbers. It may readily 
be seen that Eq.(12) is the expectation value of Eq.(13) in the coherent 
state \z) defined by a\z) = z\z) at the value \z\ = 1. This circumstance 
has been used recently to re-establish the link between the matrix element 
(z\e Xaia \z) and the properties of Stirling and Bell numbers 9 . 



2. Extending Dobinski formulas 

The purpose of this note is to show that the above results on functions of 
a^a may be extended to functions of (a^) r a s , (r, s = 1, 2, . . .) with r > s. 
Specifically, we pose the following qüestions: 

(1) What extensions of the conventional Stirling and Bell numbers occur 
in the normal ordering of [(o' ) r a s ] n ? 

(2) Can the generalised Bell numbers B r ^ s {n) so defined be represented 
by an infinite series of the type of Eq.(9) - that is, do they satisfy a 
generalised Dobinski formula ? 
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(3) May one consider the B r:S (n) as the n-th moments of a positive 
weight function W riS (x) on the positive half-axis, and may this latter 
be explicitly obtained ? 

(4) Can one attach a combinatorial significance to the sequences 
{£ r , s (n)}? 

In this note we indicate affirmative answers to the first three qüestions and 
a partial answer to the fourth. 

To this end we generalize Eq.(5) by defining for r > s: 

ns 

[(al) r aT = (a^ r ~^ £ S r ,,(n, k)(a^) k a k (14) 

k—s 

or, alternatively, 

ns 

[x r (d/dx) s ] n = x n( - r - s '> S r ,e(n, k ) xk (d/dx) k . (15) 

k—s 

Eqs.(14) and (15) intxoduce generalized Stirling numbers S rjS (n,k) which 
imply an extended definition of generalized Bell numbers: 

ns 

B r!S (n) = J2 S rAn,k). (16) 

k—s 

Note that the Stirling numbers SV,i(n, fc) have been previously studied 10 . 
Also B hl (n) = B{n) of Eq.(6). 

We have found a representation of the numbers B r ^ s (n) as an infinitc 
series, which is a generalization of the Dobinski formula Eq.(9). For r = s 
one obtains: 

n-l 

, n=l,2... (17) 
with B rjr (0) = 1 by convention. For r > s the corresponding formula is: 

oo 

B r , s (n) = [(r-sy^/e]J2 

k=0 

The formula Eq.(9) and its extensions Eqs.(17) and (18) share a common 
feature, namely, the fact that they give rise to a series of integers is by no 
means evident. 

Choosing various pairs (r, s) gives alternative representations of many 
integer sequences. Some examples are: 

(1) (r>l,*=l) 

B rA (n) = [(r - IT/e] r(n + ^)/[k\ r(^)] 



Br,r(n) = (l/e)jri 



n 

3=1 



r(i + £*) 



fc+j ' 

r— s ' 
fc+j ' 



, B r>s (0) = l (18) 
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(2) The pair (r + 1, r) leads to a hypergeometric function pF q : 



B r+hr (n) = (l/e)[f[ {U x 



3 = 1 



x r F r (n + l,n + 2, . . . , n + r; 2, 3, . . . ,r + 1; 1) 

(3) as does (2r, r) 

S 2 r,r(n) = [(m)!/er!]iFi(m + l;r + l;l). 

A still more general family of sequences arising from Eq.(18) has the 
form (p, r = 1,2,...): 



>pr (n) = {í/e) 



n 



x r íV(pn+ 1, . . . ,pn + 1 +p(r - 1); 1 + p, . . . ,1 +p + p(r - 1); 1). 
For example, for p = 3, r = 2, this reduces to 

*,«(») = (Ve) = ^ + ^ = ^ + ^ 2 F 2 (3n + 1, 3n + 4; 4, 7; 1)(19) 

whosc first four tcrms arc 1,207775,566828686621,9011375448568566265. 

Knowledge of the generalized Stirling numbers in Eq.(14) solves the 
normal ordering problem for [(o^) r a s ] n . We are able to give the appropriatc 
generating functions for the sequences B rs (n). It thcn follows that, at least 
formally, we can furnish the generating functions for S r . s (n,k) as well 11 . 
Additionally, it turns out that in certain circumstances one may obtain 
explícit expressions for them. We quote two such cases: 



and 



S r , r (n, k) — ^ ' 

p=0 



S 2 ,i(n,k) 



k-r ( l^v\ 1" 



í(k-p-r)\i 



(k — p)\p\ 



- 1 
fe! V k - 1 



(r < k < rn) 



(í<k<n) 



(20) 



(21) 



which are the so-callcd unsigncd Lah numbers 3,10 . 

For those pairs (r, s) for which we have an explícit expression for S r ^ s (n) 
we may generalize Eq.(13) to obtain the normal ordered form of e A ( at ) raS . 
For example, the matrix element (z\e x ^ ar ^ a \z) leads to the following nor- 
mally ordered expression: 



,A(a t ) r a\ 



) =: exp{[(l - Xiaïy-^r - 1))^ _ i] t } : ( 2 2) 
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We apply to Eq.(22) the same method whereby we obtained the EGF of 
Eq.(12) by taking the expectation of Eq.(13) in the coherent state \z), to 
get: 

(z\e x ^ )ra \z) =: exp{[(l - X(z*) r ~ 1 (r - 1))^t _ 1] : (23) 

which evaluates at z = 1 to give 

{z\e x ^ ra \z) x=1 =: exp{(l - A(r - 1))^ _ i} : ( 2 4) 

which is precisely 10 the EGF for the numbers B rt i(n) . 

For general r > s the corresponding S r>s (n) grow much more rapidly 
than n\ and thus may not be obtained via the usual form of EGF. One 
instead defines the EGF in terms of B riS (n) / (ni)* where t is an integer 
chosen to ensure that J2^=o B r , s {n) / {n\) t+1 has a non-zero radius of con- 
vergence. As a rcsult one obtains variants of Eq.(23) involving different 
hypergeometric functions 11 . 



3. Generalized Bell numbers as moments of positive 
measures 

The formulae Eq.(17) and (18) can be used to demonstrate what we consider 
to be the main thrust of our results, namely that B r s (n) is the n-th moment 
of a positive probability measurc. For r = s this takes the form of a sort of 
Dirac comb, while for r > s gives a continuous distribution. Consider the 
case r = s first. We rewrite Eq.(17) as 

oo 

B r , r (n) = (1/e) £ [k(k + 1) . . . (k + r - 1)]" n = 1, 2 . . . (25) 

k=o ^ >' 

which immediately indicates that B r ^ s (n) is the n-th moment of 

«m./..e»|^ (26) 

k=0 

which is a "rarefied" Dirac comb whose delta spikes are situated on the 
x-axis at x = k(k + 1) . . . (k + r — 1) fe = 0, 1, . . .. 

For r > s it is necessary to excise the part of formula Eq.(18) which 
contains n, that is Y\j =1 T(n+ for fixcd k. We reparametrize this with 

n = a — 1 and consider it as a Mellin transform; that is, we seek t^fi (x) 
such that when extended to complex a 



L 



OO * T i 

/-'^l (x)dx = J] r(<7 - 1 + —t). (27) 



3=1 



7 



The positivity of Wrfs {%) follows from the fact that it may be considered as 
resulting from an s-fold Mellin convolution of shifted Gamma-functions 12 ' 13 
which does preserve positivity. We conclude that w^s {%) is positive and so 
is W rvS (x), where 



f 

Jo 



x n W r , s {x)dx = B r ^ s (n). 



(28) 



Use of the inverse Mellin transform yiclds many analytic forms of 
W rtS {x). We exhibit some solutions of the Stieltjes moment problem 
Eq.(28): 

(1) Forr = 2,3,... 



W rA (x) = 



e(r — 1) \r — 1 



k=0 



k\ V r - 1 



exp 



/r 



X 



r-l 
r + k 



1 



(2) In particular 



W 2 ,i{x) =e~ x h {2y/x)/{ey/x), 



where I\ (y) is the modificd Bcsscl function of the first kind. 
(3) Taking r = 3 and s — 1 we have 



W 3 ,i(x) = 



1 



e\/8x 



e 2 



— \ ( 1%-) ^uF 2 (-.2:- 



T T 8 



V2 



(4) whilc for (r, s) = (5, 2) we obtain 
W 5 ,2(z) 



2^^ ( ¥ ) ·^ 2(X) ' Wh6re 



U 5 ,2Í X = 



32. 24VÏÏoF4 
+8 • 3«a;* oF 4 

+3 ■ 3^x* 0F4 



+ 
+ 



1 2 4 5 a; 
3' 3' 3' 3' 24Í 

2 4 5 6 x \ 
3' 3' 3' 3'243/ 
4 5 6 7 x V 
3' 3' 3' 3' 243/ 

and K v {y) is the modificd Bessel function of the second kind. 
(5) Finally 

1 2-3r i J_ 

^rrW = £ 2 ^ e Xr I r (2x^). 

e r 
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Our proof of the existence and positivity of W TtS (x) implies the following 
Theorem: For all the sequences defined by Eqs. (17) and (18) the following 
inequalities are satisfied: 

det [B r>a (i + j - 2)i< ÍJ <„] > (29) 
det[B r , s (i+j- l)i<i,j<„] > forall n = l,2.... (30) 

The above expresses the positivity of the two series of Hankel-Hadamard 
determinants spanned by the moment sequences. This positivity is the 
nccessary and sufficient condition for the existence of a positive measure 
W r}S (x) resulting from the given sequence of moments 14 . 

For completeness we give here the two asymptotic expansions as n — > oo 
for the series B^\(n) and B^^\(rï), obtained by the method of Hayman 15 . 
B 2 ,i(n) = (z\[(a*) 2 a] n \z) I=1 has the EGF 

and has the following asymptotics: 

B 2 i(n) " -J=(n"i + -U"! + 0(n-i)) ■ n n cxp(-n + 2y/K). (32) 
V2e 12 

Similarly, for B 3A (n) = (z\ [{a^fa] n \z) z=l 

^ = f^" (33) 

n=0 

oi i 

B 3 ,i(n) ™ — (n-3+2-?n-i+0(n- 1 ))·(2n)"exp(-n+(3/2)(2n)3).(34) 
V3e 

The role played by these approximations is analogous to that of the Stirling 
approximation to n-factorial. The growth of other B r ^ s (n) with n is in gen- 
eral even more ràpid. However, we do not have the asymptotic expansions 
for thcm at the present time. 



4. Postcript 

Our interest in integer combinatorial sequences (see also 11 ' 12 ' 17 ' 18 ) is a con- 
sequence of our parallel work on the construction of complete sets of co- 
herent states (sce 16 and references thercin). Completeness - or, equiva- 
lently, the resolution of unity property - is an essential ingredient in such 
construction, and can be shown to be satisfied if one can furnish an ap- 
propriate positive measure W(x). Thus consider a state \z) constructed 
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from an arbitrary discrete set of states \m) which are both orthonormal 
((m|m') = óm^n') and complctc: 

\z) = N~H\z\ 2 ) J2 "TT^ m >' p ^ > °" ^ 

m=0 VpM 

Resolution of unity is achieved via 



z)^(|z| 2 )(z|=/^^|m)(m|. (36) 

m=0 

The conditions Eqs.(35) and (36) boil down to a relation closely resembling 
Eq.(28), namely 

~ W{x) ~ 

_N(x)_ 

where it is assumed that the normalization of N(x) defined by 



Jo 



00 



dx = p{n), x=\z\ z n = 0,l... (37) 



has a finite radius of convergence. It is clear from Eq.(28) that when idcnti- 
fying W riS (x) with ttW(x)/N(x) we can use the B r ^ s {n) as p(n) to construct 
states of the form 



oc 



\z) r ,s = NrJ(\ Z \ 2 ) J2 /o-rr l") (39) 



where 



n—O ' • ' 

As pointed out above, B rs (n) grows much more quickly than n\ so the 
series Eq.(40) is very rapidly convergent for all r,s. This ensures that the 
problcm defined by Eq.(37) is well-defined for all x > 0. 

The properties of the coherent states defined by Eq.(39) remain to be 
investigated but they automatically satisfy the resolution of unity, due to 
Eq.(28). 

We are currently attempting to unravel the combinatorial meaning of 
some of the B r ^(n). It appears that this can be done systematically for 
B r ^{n) at least, which will be the subject of a fortheoming publication 19 . 

The integral representations of combinatorial numbers, used by us for 
the construction of coherent states, have triggered off applications for hyper- 
determinants and Selberg integrals 20 . Although of evident utility, Stirling 
and Bell numbers have not yet found their way into most textbooks on 
mathematical physics, a notable exception being that of Aldrovandi 21 . 
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